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Abstract 

The calculation of the 4 He energy and specific heat is carried out in a wide temperature 
range within the two-time temperature Green functions approach. The approximation 
improving the random phase approximation is developed providing the correct behaviour 
of the calculated specific heat at the temperatures < T/T c < 0.65 and T/T c > 1.32 
where T c stands for the phase transition point. 
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1 Introduction 

The problem of the theoretical description for the many-boson system thermodynamics in a 
wide temperature range is a very curious one and still remains unsolved. Liquid 4 He is the best 
example for the theretical studies because of its anomalous low-temperature behaviour. It is 
well-known how to obtain the thermodynamic functions in the limit T —>■ (this is given by 
the theory of Bogoliubov [l|]). The renormalization group theory also gives a fairly well results 
at the point of the lambda-transition H. But still, quite a large temperature range might be 
studied by means of numerical computations only. 

Sears || utilized neutron scattering results in order to determine the internal energy of 
liquid helium-4 via some simple thermodynamic relations. For these calculations, the potential 
of Aziz et al j|] was used for representing the He-He interatomic interaction. 

Path Integral Monte Carlo (PIMC) technique which allows for the study of the physical 
quntities temperature dependences was utilized by Ceperley and Pollock || for the helium-4 
thermodynamic properties calculation at the temperatures above 1 K. Kinetic energy of liquid 
and solid helium-4 was calculated in |J at high temperatures. Such numerical simulations are 
numerous and we do not intend to analize them here. The overview of Path Integral Monte 
Carlo (PIMC) techniques application in the theory of condensed helium is adduced in the paper 
by Ceperley 0. 

The derivation of the analytical relations for the thermodynamic properties was made by 



Vakarchuk []1T| . In the first of the cited works, the expression for the many-boson two-particle 
density matrix was written. The second work contains the analysis of the non-linear density 
fluctuations contribution into the thermodynamic properties of the Bose-system. The expres- 
sions found there are expected to provide a good description in a wide temperature range. 

In this work, we will proceed from a very well studied random phase approximation (RPA) 
which allows one to receive a low-temperature asymptotics for the thermodynamic and structure 
functions of liquid 4 He. Then, making some quite natural assumptions, we will develop what 



might be called 'effective post-RPA' since the post-RPA corrections will be taken into account 
in an effective manner. In principle, such a technique will move us substantially ahead form the 
low-temperature region. But unfortunately the way of the equations solution applied here does 
not allows to grasp the temperatures \T/T C — 1\ < 1/3 where T c stands for the lambda-transition 
point, this seems to be rather a technical problem. 

In this paper, we use the previously obtained potential [§] for the calculation of thermody- 
namic properties of liquid helium. A two-time temperature Green function technique will be 
utilized for the calculation of the energy and specific heat. The application of this technique is 
well-known, it was used by Tserkovnikov |§ for the non-ideal Bose-gas in the random phase ap- 
proximation. Additionally, this method leads to a quite short way to the RPA (Bogoliubov's) 
results. Unlike what is usualy done, the collective variables being the density fluctuations 
Fourier components are utilized for constructing the Green functions instead of the creation- 
annihilation operators. 

The organization of the paper is as follows. In Section II, we introduce the operators for the 
Green functions and write the Hamiltonian. In Section III, the equations of motion are written 
and solved. Next, the function 7fc(T) for the accounting of post-RPA terms is introduced. The 
equations for it are adduced in Section IV. Also, the expression for the free energy is written and 
other thermodynamic functions are derived from it in this Section. The discussion is presented 
in Section V. 



2 Collective variables and Hamiltonian 

We consider operators 
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the operators pk being collective variables, operators tx^ are usually known as 'currents'. No 
special name will be given to £k- 

The Hamiltonian H of the system is written as H = H^ + AH where corresponds 
to RPA and AH is the respective correction. Hereafter we will consider only RPA term and 
therefore use H to designate H^: 

H = K + $ (4) 

where K and <f> being the kinetic and potential energy in RPA respectively. 

We have the following expression for the potential energy operator in p-representation: 
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where means the Fourier transformation of the interatomic potential in helium ||, while the 
kinetic energy is given in a usual way: 

fc2 N 

^ = -^E v ?- (6) 

3=1 



It is interesting as well to consider the p-representation for K. Neglecting terms with two sums 
over the wave vector (i. e., leaving the RPA term only) we obtain 



(7) 



e k is the energy of a free particle: 
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3 Equations of motion for Green functions 

Since the total energy of the system E is the average value of the Hamiltonian we need to know 
the way of averaging our expressions. In order to find the average values of such expressions as 



PkiPk 2 ? ^kiPk 2 ) etc. we will use two-time temperature Green functions [1C 



((A(t)\B(t>))) = i9(t-t'){[A(t),B(t'))) (8) 

with operators given in the Heisenberg representation, 9 is the Heaviside step function. The 
following designations will be used: 

G pp (k 1 , k 2 ) = ((p kl | p k2 )), Gpxfa, k a ) = ((p kl 1 7Tk 2 )), • • • • (9) 

The averaging might be symbolically represented by an operator A : 

(AB)=AG BA . (10) 

We consider here only the case of k 2 = — ki = k since such expressions appear in our 
Hamiltonian ((^), (^). Therefore, the indexing of ki, k 2 will be skipped for the sake of simplicity. 

3.1 p- representation for the kinetic energy 

It is easy to show that one can obtain the following equations of motion 
hjG pp (u) = -£ k G pp {u) + G vp (u), 



huG np (uj) = £ k G np (uj)+el(al-l)Gpp(uj) + -e k , (11) 
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where 



a k = ^l + — Q = -. (12) 

Using the solutions of the above system we receive the following mean values 
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Here, the average (pkP-k) is nothing but the structure factor of the system. 



(p-kpk) 


= AG pp (k) 


= s k 


(P-k^k) 


= AG np (k) 


= n k 


(TT-kTTk) 


= AG nw (k) 


= n k 



3.2 Exact kinetic energy 



In this case, the equations of motion read: 

fuoG pp (u)) = -e k G pp (uj) + G np (u), 

hu G np (u) = -e k G np (u) + e\(a\ - 1) G pp {uj) + + -e k . (14) 

7T 

Unlike the previous case when the kinetic energy operator was given in the p-representation, 
now it turned to be imposible to solve this system with respect to G pp even in the RPA. The 
problem lies in the appearence of the functions with the higher orders of (— i/?, 2 V/m) in each 
equation following the equation for G^ p , etc., compare with (|TT|). Therefore, one should apply 
some approximation in order to obtain the closed system of the equations of motion. 

Let us consider the system of two equations for the functions G pp and G^p making a sug- 
gestion that 

G ip - 2s k G np = 2 lk (T)e 2 k G pp , (15) 

with the factor of 7&(T) being of the energy dimension. Its form will be choosed further. 

Such a substitution returs us to system flTT| ) but the quantity a k looks different now, we 
will refer it as a k . 
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One should keep in mind that a bar over a letter means 'exect with respect to RPA'. 

3.3 General properties of 7fc(T). 

Let us consider the cases of low and high temperatures. 

In the low-temperature limit the results obtained in p-representation are correct (the theory 
of Bogoliubov). Thus, 

7k(T) = 0. (17) 
In the high-temperature limit, one has the expression for the structure factor: 

S h (T) = - 1 ■ (18) 
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It leads to 



2 7fc (T) =2T = ~K, (19) 
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where K is the kinetic energy per particle, in the high-temperature limit it equals to the classical 
expression K = 3T/2. 

4 Equation for 7fc(T) and calculation of energy 



One can calculate the free energy of the system using the expression [11 



F = F + I dA$(A), (20) 
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Fo is the free energy of non-interacting system, A is the 'interaction turn-on' parameter, $ is 
the potential energy given by 
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In order to integrate $(A) in (|2T)| ) we will use the integration over a k instead of A. This 
involves the derivative d'jk/dX. One cannot found it without knowing the explicit dependence 
7fe(A). In this situation, the most simple way is to use the method from Jl2| and to write the 
so called first difference: 
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since A A = 1. Here, 7° is the quantity 7^ for the ideal Bose-gas. 
4.1 Low temperatures, T <T\ 

A self-consistent relation for obtaining the quantity 7 9 (T) might be sketched as follows: 
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One should remember that every RPA-expression contains a k with 7fc(T). 
The left-hand side of the equation (E3T) equals 
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where relation (|T^) is taken into account. 
One has for the free energy: 
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As a result, relation (E23|) is rewritten as follows: 
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where primes mean derivation over {3 and 



* = + » (27) 
For the a k , a° derivatives we obtain: 

"fc = -— > a k = —o- ( 28 ) 

The quantity 7° still remains unknown. We will try to find it utilizing the thermodynamic 
relations for the ideal Bose-gas. 

Let us consider an arbitrary (no matter, interacting or non- interacting) system first. Its 
structure factor in the long-wave limit reads: 

1/3 T 

SJT) = lim — coth — EkOik = • (29) 
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The expression e obviously equals to zero but we reserve it in the sence of the limit k —* 0, 
this will be necessary for the ideal gas case. Therefore, 

T 1 
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where the connection between structure factor at k = and the isothermic compressibility 
is taken into account. For the ideal gas V\. = 0, thus at low temperatures one receives: 

18 tor T<T " (31) 

where P$ is the ideal Bose-gas pressure and T c is the Bose-condensation temperature. As a 
result, the ideal Bose-gas structure factor S® reads: 

S° k (T) = f oc 1 (32) 



as it should be [1^ 



Now, we accept for the simplicity that the dependences on the wave vector and the temper- 
ature in quantities 7^ and 7° separate: 

7fc (T)= 7 (T)/ fc , 7 °(T) = !\T)fl f = / ° = 1. (33) 

Te above statement means the separation of the 'primary' temperature dependence in the 
comparison with the 'secondary' dependence on the wave vector. The latter will be integrated 
up when calculating the thermodynamic properties. The condition /o = /q = 1 is nothing but 
the specific normalization condition. 

Now, from 7$ = one obtains 7 = and, taking ( p%p into account, also 7°' = 0, a®' = 0. 
These equalities are valid for T < T c . 

Let us consider the expansion of ( fffi ) in terms of 7 in order to obtain the equation for 7fc(T). 
The first-order terms read: 
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where the dot denotes the derivative with respect to the temperature and a q is defined by the 
Bogoliubov's expression (|T2| ) 

Similar expansion for energy (123) leads to the following equation: 



E 



4tt 2 q 



dk k 2 e k 



a k coth 



2T 



a 2 k + l 



-1 



Att 2 q 



+ Y 



An 2 g 



dk k 2 f , 



'i— 

\ATa k 



+ 



e k coth ^ 



2T 



T 



sinh 



2 aoi 



dkk 2 f k —e k 



+ 



1 
1 



2T 



2Te k a k " 



sinh 



coth 



2 e k a 



k^k 
2T 

2T 



(35) 



The last item in the second integral causes some difficulties since it does not converge at 
k — > oo. One can ensure the convergence by means of the function f k (it is arbitrary by 
now). The next problem with this item is the incorrect behaviour of the heat capacity at low 
temperatures. As it will be shown further, j(T) is proportional to T when T —>■ 0. It leads to 
the linear temperature dependence of heat capacity while one should have oc T 3 . We can avoid 
this problem by choosing in the special form. We will demand the following: 

1 
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The similar problem is with giving quadratic dependence with respect ot T. 
second equation: 
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Only a complex f k can satisfy both these equation. For the sake of simplicity in calculations 
we have choosed f k in the following form: 

< k < ki, 
k\ < k < &2, 

k 2 <k< k 3 , (38) 
k% < k < k^, 
otherwise 
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with equally spaced intervals, k n = n x 2 A -1 . Three quantities fn, fm, fw are necessary to 
satisfy both real and complex terms in the equations (|36|-p7|). 
The equation for 7 might be rewritten as follows: 

a 1 (T)T 7 + a (T) 7 + O(7 2 ) 
Neglecting higher-order terms one obtains: 
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where the constant C is unknown yet. One can define it studying the behaviour of the isothermic 
compressibility near the absolute zero (30). We received the following value by averaging the 
results of the data Jl4| interpolation at the lowest temperatures over four, five and six points: 

C = 7 (0)~6.6. (41) 

The j(T) curve is adduced at Fig. [I]. 



One should notice that teh correct results must be expected only for the temperatures at 
which it is possible to neglect the higher-order terms in the equation for 7 ^ (|34]) . In other 
words, if 2^ k /e k > 1 then the results are not correct. Taking into account that the typical 
wave vector absolute value for 4 He is 2 A (it corresponds to the first maximum on the structure 
factor curve) we obtain the upper limit for the applicability of our results as T/T c ~ 2/3. 




Figure 1: The dependence of 7 on the temperature. 



4.2 High temperatures, T > T\ 



In this case one can use the techniques similar to the described in the previous sections for 
low temperatures. Unfortunately, this seems to be a rather complicated problem. Firstly, the 
function 7°(T) is not equal to zero already. Secondly, the quantity jh(T) is not small in this 
temperature region and it is not possible to use the expansion into series, we also did not 
manage to solve equation (^) (the initial condition for the differential equation is unknown). 
Therefore, for the simplification of the calculations we decided to substitute ^k{T) with 2K/3, 
see ([TDD, where K is the kinetic energy per particle. 

We have rewritten the expression for the free energy ( p5| ) in order to obtain the equation 
for AK = K — K : 
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where the dot means derivation with respect to T: K = dK/dT. 

Having calculated the kinetic energy form Eq. ( [42]) we then proceed to the calculation of 
the total energy. After the numerical derivation of the latter we receive the specific heat of the 
liquid helium. This result is adduced in Fig. [| together with the data for the low temperatures 
from the previous subsection. 



Figure 2: Energy of liquid helium-4. Solid line — calculated energy, dotted line shows the 
regions where the approximations used cannot be accepted; dashed line — energy of ideal 
Bose-gas; dashed-dotted line — Bogoliubov's approximation. 



5 Discussion 



As one can see from Fig. ^| we obtained a fairly good description in the region < T/T c < 0.65, 
as it was expected. In the normal phase (T > T c ) a similar situation is received for the 
temperatures T/T c > 1.32. It is necessary to say that the phase transition temperature T c 
in our case coincides with the respective temperature of the ideal Bose-gas. One should use 
the effective mass m* of the helium atom instead of that of free atom m = 4.0026 a. m. u. 



As will be shown elsewhere fL5 |, it is possible to receive the effective mass in a sel-consistent 



manner within the Green function formalism. The ground-state (T = K) value equals to 
m* ~ 1.58 m. We assume that the temperature dependence of the effective mass is significant 
in the temperature region under consideration: m*(T) ~ const = m*(0) = 1.58m. Therefore, 
one receives the transition temperature of 1.99 K vs experimental 2.17 K ||15||. 
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